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Analytical all-orders results are presented for the one-renormalon-chain contributions to 
the Bjorken unpolarized sum rule for the F\ structure function of uN deep-inelastic scat- 
tering in the large- Np limit. The feasibility of estimating higher order perturbative QCD 
corrections, by the process of naive nonabelianization (NNA), is studied, in anticipation of 
measurement of this sum rule at a Neutrino Factory. A comparison is made with similar 
estimates obtained for the Bjorken polarized sum rule. Application of the NNA procedure 
to correlators of quark vector and scalar currents, in the euclidean region, is compared 

■ ! with recent analytical results for the 0(ajN F ) terms. 
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1 Introduction 



The differential cross-section for neutrino (anti-neutrino) nucleon deep-inelastic scattering 
(DIS) is parametrized by three structure functions in the familiar formula 
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with < x < 1, y — Eh a d/E u , < y < 1/(1 + xM-^/2E u ), and a +(— ) sign applying to 
neutrino (anti-neutrino) beams. 

Analyses of data from previous uN DIS experiments concentrated on the extraction of 
F 2 and xF 3 (see e.g. the works of Ref.gl]) and the ratio R = (1 + 4M^r/Q 2 )F 2 /(2:rF 1 ). 
(For a recent extraction in uN DIS, see Ref.[[|.) However, there are already at least two 
attempts to obtain direct information about F± N ||. Such efforts are interesting from 
several points of view. They allow comparison with models for the 0(1/Q 2 ) contributions 
to both F\ and xF%, obtained within the framework of an infrared renormalon approach, 
which predicts similar x-dependence of these two power corrections Q. A second point 
of interest is comparison with QCD sum rule estimates for the twist-4 matrix element 
introduced in Ref.|| which gives an 0(1/Q 2 ) correction to the unpolarized Bjorken sum 
rule (Bjunp SR) 
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with a corresponding twist-4 term for the Bjorken polarized sum rule estimated in Ref.j7|. 

In this note, we present analytical all-orders results for the one-renormalon-chain con- 
tributions to Cejunp, obtained within the large- Np expansion, and then apply the simplistic 
procedure of naive nonabelianization (NNA), proposed in Ref. [|J. A comparison is made 
with the polarized Bjorken sum rule, and with the application of NNA to vector and 
scalar correlators. We hope that these considerations will encourage theoretical study of 
the Bjorken unpolarized sum rule, whose experimental study may be enabled by uN DIS 
data from a future Neutrino Factory, operating in a region of medium Q 2 , with iV^ = 3, 4 
active flavours. (For active consideration of such a facility, see Ref. 0.) 



2 Large- Np series 



Here we study radiative corrections to Bjunp SR induced in the large- Np limit. Using 
methods developed for the calculations in Refs. |l~0|, [TT| , [i~2| , the large- Np limit of the 
perturbative part of Eq.(2) was obtained with the result: 
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where a s =a s (fi 2 = Q 2 )/An is the coupling in the MS-scheme, Cf = 4/3, Tp = 1/2 and 
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This expression produces the following large- N F MS-scheme series for the Bjunp SR 
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531441 4782969 43046721 v ; 

It should be noted that in this new result the order x 2 and x 3 -terms are in agreement with 
the corresponding parts of the total expression for the a^-correction to the Bjunp SR |13| 
and its a^-term, calculated in the MS-scheme in Ref.[Q. Notice also that the series of 
Eq. (5) has sign-alternating behaviour with factorially increasing coefficients. This pattern 
is explained by the dominant role of the renormalon at 5 = 1, which is generated by a 
single chain of quark-loop insertions into the corresponding one-loop QCD diagrams for 
this characteristic of uN DIS. 

It is interesting to compare the above expressions to the analogous ones, obtained in 
Ref. fllT] , for the coefficient function C*Bj P (Q 2 ) in the Bjorken polarized sum rule (Bjp SR) 
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The large iVp-expression for the Bjp SR can be obtained from the following equation JIT 
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The all-orders large- Np result of Eq.(8) was given in Ref.pi|. The expansion up to order 
0(a®Np) was also given in Ref. |lT| 



in the MS-scheme and reads: 
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As in the case of Eq.(5), this series is dominated by the renormalon at 5 = 1, which has the 
same reside in each sum rule. (Note that the renormalon in K(S) at 5 = —1 is suppressed 
by a factor 7jexp( — 10/3) = 0.018, relative to the dominant renormalon at 5 = 1.) 

In the next section we apply naive nonabelianization to the large- iVp series for both 
unpolarized and polarized Bjorken sum rules, in an attempt to estimate higher-order 
perturbative coefficients. 
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3 Application of NNA to euclidean sum rules 



By naive nonabelianization (NNA) we simply mean the substitution || Np — > — |/3o = 
Np — y in the leading terms of the large- Np expansion. The hope is that terms of lower 
order in Np may be roughly estimated by assuming that they follow the leading terms with 
weights generated by the one-loop term of the QCD beta function, /3q = 4- Ca — ^TpNp. In 
a variety of cases (see e.g. Refs.||, 0, |l|, [L6| |17| ) this simplistic procedure gives reasonable 
estimates of known higher-order perturbative coefficients in different physical quantities 
(see Ref -11181 for a review). As a rule, the signs of contributions from sets of diagrams with 



fewer quark loops are correctly predicted and often the actual magnitudes of coefficients 
of lower powers of Np are within a factor of 2 of the NNA estimates. In this section we 
apply this procedure to the Bjorken unpolarized and polarized sum rules, where Q 2 > 
is, of course, in the euclidean region. 

We expand the perturbative contributions to Eq.(2) in powers of a s /n = Aa s 

CBjunp = l + E^f-) n ( 10 ) 
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with known results d\ = — g and 



d 2 = -3.8333 + 0.2963(W F 

d 3 = -36.155 + 6.3313iV F - 0.15947iV| (11) 



from Refs.[|T^, II], in the MS-scheme. Taking the input of Eq.(5), naive nonabelianization 
yields 

d NNA = _4.8889 + 0.29630A^ F 

^nna = _43_4 14 + 5.2623iV F - 0.15947iY| 

^nna = _457.02 + 83.094A^ F - 5.0360A^| + 0.10174AT| (12) 

obtained from the exact results for the terms with highest powers of Np. One can observe 
reasonable agreement between the estimates and exact results for the coefficients of N F in 
g?2 and of Np in d%. Moreover, even the estimate, —43.414, of the Np term in d 3 has the 
correct sign and a magnitude close to the known value, —36.155. 

We now perform a similar analysis for the Bjp SR perturbative series Cb] P (Q 2 ) = 
1 + J2n>i d n {ot s /ir) n . The exact known results are di = —1 and 

d 2 = -4.5833 + 0.33333iV F 

d 3 = -41.440 + 7.6073^-0.17747^, (13) 



obtained at 0(a%) in Ref.pH and at 0(a 3 s ) in Ref.pj. The NNA estimates are 

—NNA 

d 2 = -5.5 + 0.33333A F 

—NNA o 

4 = -48.316 + 5.8565iV F - 0.17747^ 

^nna = _ 466 _oo + 84.728A^ F - 5.1350A^| + 0.10374A^| (14) 
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with g?2 and in tolerable agreement with the exact results of Eq.(13). 

More ambitious estimates are presented in Tables 1 and 2, where we have evaluated 
Eq.(12) and Eq.(14) for Np = 3,4,5 active flavours and compared these estimates with 
exact results from Eq.(ll) and Eq.(13), or with estimates in Ref.pl| at O(a^), obtained by 
a method of effective charges |22|. (The latter are in fair agreement with Pade estimates 
in Ref. @.) 
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Table 1. The A^-dependence of the NNA expressions for Bjunp SR and their comparison 
with the results of explicit calculations and 0(ot\) estimates of Ref . |f2"T[l . 
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Table 2. The Af^-dependence of the NNA expressions for Bjp SR and their comparison 
with the results of explicit calculations and 0(a\) estimates of Ref. |2l"H . 



Here one is pushing NNA rather hard, since there are substantial cancellations between 
powers of Np with alternating signs, which are most pronounced at Np = 5. Nevertheless, 
there is agreement to within a factor of 2 for Np = 3, 4 at order a 2 s and a 3 . 

The similarities between Tables 1 and 2 are rather striking: the NNA estimates, the 
known values, and the effective-charges estimates in the two sum rules follow very similar 
patterns. Indeed this was observed in the renormalizat ion-group invariant analysis of 
known results in Ref.pl], where this similarity between polarized and unpolarized sum 



rules appeared to be somewhat mysterious. From our point of view there is a simple- 
minded explanation: the residues of the dominant infrared renormalon, at 5 = 1, in the 
Borel transforms of Eq.(4) and Eq.(8) are identical. Thus one may attribute the close 
similarities in the full perturbative structures to the rough success of NNA, which assumes 
that the overall trends are driven by this renormalon. 

However, in the concluding section we remind the reader of two important effects that 
suggest caution in relying on estimates obtained from a single chain of quark loops. 
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4 Application of NNA to euclidean correlators 



To sharpen our understanding of the limitations of NNA, we reconsider analyses of the 
vector and scalar ]TJ| correlators, of light-quark currents, following the recent and 



impressive calculations of the 0(a^N F ) terms reported in Ref. |25 
In the vector channel we study the Adler function 
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where R(s) is the well-known e + e ratio, Qf are the quarks charges, and the known 
perturbative results, at large Q 2 > in the euclidean region, are dj = 1 and 
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taking d% from Ref.p6f and d^ from Ref. ||27||, with neglect of the term in the vector 
correlator with two quark loops and a three-gluon intermediate state. The 0(a A s Np) term 
in dj was obtained in Ref.|28[, and the 0(a A s Np) term was recently published in Ref. f25fl . 
The corresponding NNA estimates are 
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While the pattern of agreement of signs and rough magnitudes is comparable to that for 
the sum rules it is notable that the 0(Np) term in d% exceeds the NNA estimate by a factor 
1.5 and the O(Np) term in dj exceeds the estimate by a factor 1.8753/0.49959 ~ 3.8. 

Proceeding to the scalar correlator, again in the euclidean region, we study the following 
analogue of Eq.(15) 
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648.71 - 63.742iV F + 0.92913iV| 

d lo + dt^Np + 54.783iV| - 0.45374iV| 



(19) 



taking <if from Ref. |]29[| , <if from Ref. |30[| , the 0(a^Np) term in d% from Ref. Hl2| , and the 
0(a i s Np) term from Ref.[pj|. We choose this comparator so as to make contact with the 



analyses of Refs.[31, EH], notwithstanding the comment made in Ref. [12] that the necessity 
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of a second subtraction in the dispersion relation for the scalar correlator makes such a 
construct infrared unsafe. The corresponding NNA estimates are 



d s,NNA = 3i.4 6 5 _ i.9070iV F 

^s.nna = 252.96 - 30.661iV F + 0.92913A^| 

^s.nna = 2038.3 - 370.60iV> + 22.460iV| - 0.45374iV| (20) 

again with a fair pattern of agreement in signs and magnitudes. The O(Np) term in df 
exceeds the NNA estimate by a factor 2.1 and the 0(N F ) term in c?| exceeds the estimate 
by a factor of 54.783/22.460 « 2.4 

Two remarks are in order. First, one notes that the absolute sizes of the radiative 
corrections in the scalar channel are much larger than those in the vector channel. This 
may be regarded as a success for the NNA estimator, which attributes this trend to the 
inadequately subtracted dispersion relation of Eq.(18), taken as the object of study in 
Refs. f3"I|, f25| . The failure to make a second subtraction in the scalar channel results |12| 



in an infrared renormalon at S = 1, whereas no such renormalon can appear in Adler's 
correctly subtracted dispersion relation of Eq.(15). Moreover, this disparity between the 
scalar and vector channels is already apparent in the 0(a s ) diagrams, into which quark 
loops are inserted, where the coefficient df — is almost 6 times larger than d\ = 1. 

As in the previous section, this renormalon analysis has value only if the known radia- 
tive corrections follow the pattern suggested by NNA. Fortunately this is again the case, 
in the euclidean region. We remark that the very large 0(a^N F ) coefficient d^ 2 = 54.783 
is better approximated by NNA than is the far smaller coefficient dj 2 — 1-8753. It thus 
appears that NNA is a useful indicator of significant euclidean radiative corrections. 

The second important observation is that a very different picture emerges if one chooses 
the comparator R s (s), which is the discontinuity of the scalar correlator across the cut 
in the minkowskian (timelike) region — Q 2 = s > 0. The modest success of NNA in 
the euclidean region does not ensure comparable success in the minkowskian region, be- 
cause of the numerically large terms involving powers of 7r 2 and coefficients of the beta 



function, (3(a s ), and the mass anomalous dimension, 7 m (a s ), noted in These easily 
computed effects of analytic continuation do not naively nonabelianize. Moreover, the 
infrared renormalon at 5 = 1 in Eq.(18) is absent [|l^] from the scalar imaginary part 
R s (s). Accordingly, NNA behaves poorly in this comparator. For example, it was noted 
in Ref.|25| that the exact 0(a A s N F ) term in R s (s) exceeds that obtained by applying NNA 
to the exact 0(a 4 s N F ) term of Ref.[12|] by a factor of 9.6848/1.0128 rj 10. 



Thus we learn from the analytical calculations of Refs. [12, [25| that NNA estimates for 
the scalar correlator perform better in the euclidean region, where we have already noted 
tolerable agreement with exact results on sum rules. 



5 Conclusions 



In summary: 
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1. The new all-orders results for the leading terms in the large- Np expansion of the 
perturbative contributions to the Bjorken unpolarized sum rule, when naively non- 
abelianized, offer an explanation of the observation of Ref . [[24]] that the actual radia- 
tive corrections in this sum rule closely follow those of the Bjorken polarized sum 
rule. Within the simplistic, but here relatively successful, framework of NNA, we 
attribute this parallelism to the equality of the residues of the dominant infrared 
renormalon, at 5 = 1, in the Borel transforms of Eq.(4) and Eq.(8). 

2. The pattern of NNA estimates for the vector and scalar correlators, in the euclidean 
region, also suggests a leading role for infrared renormalons. In this case, it helps 
one to understand why the actual scalar radiative corrections are so much larger 
than those in the vector case. To attribute this to NNA, it is necessary for the 
procedure to have reasonable success in each channel. In fact we find that it is 
slightly more successful in the scalar case, notwithstanding conspicuous failure after 
analytic continuation to the minkowskian region. 

3. Since NNA yields correct signs and sensible magnitudes for all known coefficients 
of terms with lower powers of Np, in all four euclidean analyses considered in this 
work, we believe that its sign predictions in Eqs. (12,14,17,20), for unknown 0{pq) 
terms, are probably reliable. In particular, we expect that in Eqs. (16,19) further 
dedicated calculation will yield negative signs for d 4 i and positive signs for the dj^ 
coefficients that result from purely gluonic radiative corrections. 

4. It is important to recognize the intrinsic limitations of NNA. First we remark that 
the Np dependence of QCD radiative corrections results not only from insertions 
of quark loops in gluon propagators; in addition there are diagrams with quark 
loops inserted at gluonic vertices. It seems unreasonable to expect NNA to mimic 
contributions such as that in Fig. 1(a) of Ref.[f2^], where a quark loop modifies a 



three-gluon vertex within a second quark loop. More significantly, perhaps, one first 
encounters at 0(af) the effects of quark loops in different gluon propagators. As 
remarked in Ref.[^2|, multiple renormalon chains modify the asymptotic behaviour 
at next-to-leading order in 1/Np, notably in the 0(aj,Np) terms. Therefore, it would 
appear more prudent to use the final row in each of the Tables 1 and 2 as estimates 
for DIS sum rules. Moreover, in the vector correlator, the effective-charges analysis 



of Ref. [pT|j appears to be in good agreement with the partial analytical results of 
Ref.ra. 
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